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barotropic Stokes systems 
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Abstract. We propose finite element methods for compressible barotropic 
Stokes systems. We state convergence results for these methods and outline 
their proofs. The principal tools of the proofs are higher integrability estimates 
for the discrete density, equations for the discrete effective viscous flux, and 
renormalized formulations of the numerical method for the density equation. 



1. Introduction 

In this contribution we consider mixed type systems of the form 

(1.1) dtg + dW{gu) ^0, in (0, T) x 17, 

(1.2) "fiAu- XDdivu + Dp{g) ^ f, in(0,T)xf7, 
with initial data 

(1.3) g\t=o = go, on fl. 

Here 17 is a simply connected, bounded, open, polygonal domain in M.^ [N = 2,3), 
with Lipschitz boundary dfl, and T > is a final time. The unknowns are the 
density g — g{t, a;) > and the velocity u = u {t,x) e R^, with t e (0,r) and 
a; G ri. We denote by div and D the usual spatial divergence and gradient operators 
and by A the Laplace operator. 

The pressure p{g) is governed by the equation of state p{g) = ag'^, a > 
(Boyle's law). Typical values of 7 range from a maximum of | for monoatomic gases, 
through ^ for diatomic gases including air, to lower values close to 1 for polyatomic 
gases at high temperatures. We will assume that 7 > 1. Furthermore, the viscosity 
coefficients /i, A are assumed to be constant and to satisfy /i > 0, A^A + 2^ > 0. 

At the boundary dfl, the system (|1.1I) - (|1.2I) is supplemented either with the 
homogenous Dirichlet condition 

(1.4) M = 0, on{0,T)xdn, 
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or with the Navier-slip condition 



(1.5) 



M • = 0, curltt X ly = 0, on (0,T) x dQ. 



System (|l.l|) - (|1.2p can be motivated in several ways. Firstly, it can be used as 
a model equation for the barotropic compressible Navier-Stokes equations. This is 
a reasonable approximation for strongly viscous fluids for which convection can be 
neglected. Secondly, in [Tj Section 5.2, Remark 5.8], Lions construct solutions to 
the barotropic compressible Navier-Stokes equations using solutions of the system 
(Hull-dill). Finally, by setting 7 = 1, / = 0, and ^ = 0, the system (IlII])-(IIl2l) is 
exactly on the same form as the model derived in 8J for the dynamics of vortices 
in Ginzburg-Landau theories in superconductivity. 

Among many others, the semi-stationary system (|l.ll) - (|1.3p has been studied 
by Lions in jTl Section 8.2] where he proves the existence of weak solutions and 
some higher regularity results. 

The plan of this contribution is to summarize some results m m [6] from an 
ongoing project to develop convergent numerical methods for multi-dimensional 
compressible viscous flow models. We construct numerical methods that comply 
with the mathematical framework developed for the compressible Navier-Stokes 
equations by Lions f7\ and Feireisl f3l. Over the years, several numerical methods 
appropriate for compressible viscous gas flow have been proposed. Except for some 
one-dimensional situations (cf. Zhao and HofF [9l llOj ). it is not known, however, 
that these methods converge to a weak solution as the discretization parameters 
tend to zero. Convergence analysis for the compressible Navier-Stokes system is 
made difficult by the non-linearities in the convection and pressure terms and 
their interaction. As a first step towards establishing convergence of numerical 
methods for the full system, we consider simplified systems that contain some of 
the difficulties but not all. In that respect l|l.ip - (|1.2p provides an example. 

The finite element methods presented here are designed to satisfy the properties 
needed to apply the weak convergence techniques used in the global existence the- 
ory for the compressible Navier-Stokes equations. Although the simplified system 
(|l.ip - (|1.2p contain additional structures rendering the solutions more regular than 
those of the full Navier-Stokes system, we strive to employ techniques that can po- 
tentially be extended to the full system. More specifically, our finite element meth- 
ods are designed such that Hodge decompositions of the velocity, u = curl^ -I- Dz, 
can be achieved and described at the discrete level. This is important since then 
a discrete equation for the effective viscous flux, (A + /i) divu. — p{g), can easily 
be extracted from the numerical scheme. It is the properties of this quantity that 
leads to strong convergence of the numerical density function; the major obstacle 
to proving convergence of a numerical method. Formally, multiplying the equation 
(|1.2p with M, integrating by parts, and using the continuity equation multiplied 
with —zjP'{q) one obtains the energy relation 



A similar relation holds for our finite element methods, which reveals the rather 
weak a priori estimates that are available to us. Indeed, it is now clear that a major 
obstacle is to obtain enough compactness on the numerical density to conclude 
that p{Qh) pis)', of course, this is equivalent to Qh — >■ g almost everywhere. 
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The remaining part of this contribution is organized as foUows: We conect 
some preUminary material, including the notion of weak solutions, in Section [2l 
In Section [3] we present a finite element method for the semi-stationary Stokes 
system in primitive variables. We state a convergence result for this method and 
comment on its proof. This method is fully developed and analyzed in [5]. In 
Section m we present and analyze an alternative finite element method [4] for the 
same system. This method is, however, restricted to the case of the Navier-slip 
boundary condition (jl.5p . Finally, we conclude this contribution by presenting a 
convergent finite element method for the Stokes approximation equations, which 
generalizes the system (|l.ip - (ll.2p by adding an additional time derivative term dtu 
to the equation for the velocity. 

2. Preliminary material 

Throughout the text we make frequent use of the divergence and curl operators 
and denote these by div and curl, respectively. In the 2D case we denote both the 
rotation operator taking scalars into vectors and the curl operator taking vectors 
into scalars by curl. We make use of the spaces 

p^curi,2(f^^ = {t; e L^{n) : cmlv £ L^{n)} , 

where v denotes the unit outward pointing normal vector on 917. If d G VF'^'^'^(ri) 
satisfies v ■ i^ldn = 0, we write v G Wj^'^'^(O). Similarly, v G VFQ"'''^(ri) means 
V G VF'^'^'^(ri) and v x i^ldn — 0. In two dimensions, xo is a scalar function and 
the space Wg"'''^(17) is to be understood as Wq''^{^}). To define weak solutions, we 
shall use the space 

W= {v e L^{VL) : divv G L'^{n),cnT\v G L'^{n),v ■ iy\on = O} , 

which coincides with W^'^'^in) n W'^'^^^in). The space W is equipped with the 
norm = Il«lli2(s:^) + ||div + WcmlvWl^^^y It is known that |H|^ is 

equivalent to the norm on the space {v G H^{il) : v ■ i^ldn — O}. 
Next we introduce the notion of weak solutions. 

Definition 2.1 (Weak solutions). A pair {g, u) of functions constitutes a weak 
solution of the semi-stationary compressible Stokes system (jl.ip - p.2p with initial 
data (|1.3p provided that: 

(1) (e, u) G i°°(0, T; L^n)) X L2(0, T; W^), 

(2) dtg + div(eM) = in the weak sense, i.e, V(/) G C°°([0,T) x Q), 

(2.1) [ [ g{(t)t + uD(i)) dxdt+ [ go(f>\t=o dx ^ 0- 
Jo Jn Jn 

(3) -fiAu-XD div u + Dp{g) = / in the weak sense, i.e, G C°°([0, T) xl7) 
for which • = on (0, T) x dn, 

(2.2) / / ji cmlu curl 4> + [{n + X) div u — p{g)] div (p dxdt = / / f4> dxdt, 
Jo Jn Jo Jn 

Whenever the Dirichlet boundary condition (jl.4p is part of the problem, we 

require that m x ^ = on (0, T) x dil in (1) and moreover that (|2.2p holds for test 

functions satisfying = on (0, T) x d^l. 
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3. A non conforming finite element method 



Following [5 , in this section we present a finite element method for the system 
(|l.ip - (|1.2p appropriate for both the Dirichlet boundary condition (|1.4p and the 
Navier-slip boundary condition (|1.5p . 

For discretization of the velocity we will use the Crouzeix-Raviart element 
space. Consequently, the finite element method is non-conforming in the sense that 
the velocity approximation space is not a subspace of the corresponding continuous 
space, W^''^{fl). Moreover, we will use a non-standard finite element formulation. 
More precisely, the formulation implicitly use the identity 



valid for all u G yV(il) satisfying any of the two boundary conditions p.Sp and (|1.4p . 
However, as the method is non-conforming, this identity does not hold discretely 
(as a sum over elements). Still, at the discrete level, the form on the right-hand side 
of p.ip is used. In contrast to the standard situation where the form on the left- 
hand side of (|3.ip is used, this discretization does not converge unless additional 
terms controlling the discontinuities of the velocity are added [T] : 



where F/i is the set of faces and |-|r denotes the jump over the edge F. 

The advantage with this formulation is that it enables Hodge decompositions 
of the numerical velocity field. By writing u = curl^ + Dz the Laplace operator 
can be split into a curl part and a divergence part plus certain jump terms. This 
is very convenient in the convergence analysis of the method. A discrete equation 
for the effective viscous flux can then be easily obtained. The reader is encouraged 
to consult 5_ for the details 

Given a time step Ai > 0, we discretize the time interval [0,T] in terms of the 
points t"^ = mAt, m = 0, . . . , M, where we assume that MAt = T. Regarding the 
spatial discretization, we let {-E/ilh be a shape regular family of tetrahedral meshes 
of f2, where h is the maximal diameter. It will be a standing assumption that h and 
At are related such that At — ch, for some constant c. For each h, let F/j denote 
the set of faces in Eh- 

We need to introduce some additional notation for discontinuous Galerkin 
schemes. Concerning the boundary dE of an element E, wc write /+ for the 
trace of the function / achieved from within the element E and /_ for the trace 
of / achieved from outside E. Concerning an edge F that is shared between two 
elements E^ and £'+, we will write /+ for the trace of / achieved from within E^ 
and /_ for the trace of / achieved from within Here and E^ are defined 
such that v points from E- to E^^ where v is fixed (throughout) as one of the 
two possible normal components on each edge F throughout the discretization. We 
also write |/]]r = /+ — /- for the jump of / across the edge F, while forward 
time-differencing of / is denoted by = - /™ and = 

We will approximate the density in the space of piecewise constants on Eh and 
we denote this space by Qh{^)- For approximation of the velocity we will use the 



(3.1) 




Wi I'^-'^M^- + I« X i^M^ X lyjr dS{x), 
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Crouzeix-Raviart [2] element space 

Vh{n) = !^VH;vh\E e V^{E), yE e Eh, ^KJ dS{x) = 0, vr e r, 

To incorporate boundary conditions, we let degrees of freedom of V/i(ri) vanish at 
the boundary. That is, for Navier boundary condition (|1.5p we require 



J vh-v ds{x) - 0, vr e r,, n dn, yv,, e v,,, 

and for the Dirichlct boundary condition (|1.4D . 

To the space Vh{fl) we associate the semi-norm 

\'"h\v^(n) = l|curl,i7;/i||^2(o) + \\ divhVh\\l2(^^-j 

+ Wr\Y. IIK-^llli^(r) + IIKx^llli2(r), 

and the corresponding norm 

\\'"h\\v,^^n) = hhWhf^n) + \'"h\v,An)- 
Here, curl/i and div^ denotes the curl and divergence operators, respectively, taken 
inside each element. The scaling parameter e > is required to prove convergence 
of the finite clement method. The size of e will affect the accuracy of the method 
and it should therefore be fixed very small in practical computations [5 . 

Before stating the finite element method, we recall from [5j the following basic 
compactness result for approximations in Vii{Q). 

Lemma 3.1. There exists a constant C > 0, depending only on the shape regu- 
larity of Eh and the size of il, such that for any f G 

and \\vh\\L^n) < C\vh\v^(n), e Vh{fl). 

Definition 3.2 (Finite element method). Let {Qhi^)} h>o ^ sequence (of 
piecewise constant functions) in Qh{^) that satisfies > for each fixed h > and 
Qh ^ ^-^^ ^ ^^'^ L^{Q) sls h 0. Set := H^/, where it is understood 
that n^/ projects f{t, x) onto constants both in time t and space x; for notational 
convenience we set /™ := fh{t^, •) & Qh{^) for any m — 0, . . . , M. 

Now, determine functions (g™, it™) G Qh{^) x Vhi^), m = 1, . . . , A/, such that 
for all (l)h £ Qh{^), 
(3.2) 

/ d1[g^;:]<t>H dx-AtY, I (^,™(m™ . v)+ + e'J^(Mr • ^r) IMr dS{x) = 0. 

and for all Vh G Vh{^), 

/ /X curlft M™ curl,i Vh + [{n + A) div^ m" - p{q^)] div^ Vh dx 

(3 3) 

+ J2tFi[ ^< ■ '^H^^. • ^^1 + K X i^Evh X i^i dS{x) = / fj:^Vh dx, 
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In (|3.2p . {uh ■ ~ max /p ■ v dS{x), o| and 

(uh ■ ly)^ (x) = min Jp it/i • dS{x), o| for x £ F and all F G F/j. 

The existence of a solution to the discrete equations p.2p - (|3.3p is proved in 
[5] by using a topological degree argument. In [5j it is also shown that the scheme 
preserves the total mass and that the density remains strictly positive provided 
that the initial density is strictly positive. Moreover, for any m — 1, . . . , M, 

m 

P{bI) dx + cY^M\\ul\\l^^^^ 

m „ m „ 

+ E / P"i0ttM-'r dx + Y.11^^ P"[g''^Mll\ui.y\ dx 
k=i-'^ fc^ireEh ''^ 




where P{g) = sMl if > 1 and P{g) = g\ogg if 7 = 1. Moreover, g'^^ e [gij; \ g'^] 
and e [(?+,£»'!]. 

Next, for each fixed h > 0, we extend the numerical solution {(fi™, it™)}^^Q 
to the whole of (0, T) x 17 by setting 

(3.4) = teitm-i,tm), m = l,...,M. 

In addition, we set £'/t(0) = g^. 

The main result of [5] is that the approximate solutions (|3.4p converge to a 
weak solution of the semi-stationary Stokes system (|l.ip - (|1.2p . 

Theorem 3.3. Suppose f e i^((0,T) x V,) and go e L''{n), if j > 1, and 
gologgo € L^{fl), if ^ = 1. Let Mft)};j>o ^ sequence of numerical solutions 
constructed according to p.4p and Definition \3.S\ Then, passing if necessary to a 
subsequence as h ^ 0, Uh ^ u in L^(0,T; L^(fl)), ghUh — ^ gu in the sense of 
distributions on (0,T) x fl, and gh ^ g o,.e. in (0,T) x fl, where the limit pair 
{g,u) is a weak solution as stated in Definition \2.1\ 

Comments on the proof of Theorem \ 3.!A In proving convergence to a weak 
solution of the continuity equation p.ip the main step is to obtain convergence of 
the product guUh gu in the sense of distributions; this follows from an Aubin- 
Lions argument using the spatial compactness of the velocity established in Lemma 
Ocombined with the fact that d^[gh] Gb L^(0, T; 

To conclude convergence to a weak solution of the velocity equation (|1.2p . we 
need a higher integrability estimate for the numerical density. We achieve this 
by utilizing test functions Vh £ Vuiyt) satisfying dvvvh ~ p{gh), thereby obtain- 
ing p{gh) Gb i^(0, T; L^(r2)). Next, we establish strong convergence of the den- 
sity. This is obtained by first proving weak sequential continuity of the effective 
viscous flux. That is, first we establish that lim^_^o [{^ + M)divith — pigh)] gh = 
{pi + A) divit — p{g)g, where the overbar denotes the weak limit. In this step, the 
div-curl structure of the scheme is utilized. In particular, we employ test functions 
Vh G Vh{^) that satisfies divvt = gh and cnvlvh = on elements away from the 
boundary. Finally, using this and a renormalized version of the continuity scheme 
(|3.2p . we obtain strong convergence of the density. 
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4. A mixed finite element method 

Following ^4], we present an alternative finite element method appropriate for 
the Navier-slip boundary condition p.5p . The method is derived by introducing 
the vorticity w — curl u as an auxiliary variable and recasting (|1.2p as 

/i curl w — (A + /i)Z? div u + Dp{g) = /, 

where also the identity —A = curl curl —_D div is used. This leads naturally to the 
following mixed formulation: Determine functions 

{w,u) e L'{0,T;Wr'-\n)) X L'{0,T;Wt-\n)) 

such that 



(4.1) 



/ / ^ curl wv + [{^ + X) diy u — p{g)] diY V dxdt — / / fv dxdt, 
Jo Jn Jo Jn 



WT] — cuArju dxdt = 0, 

o 

for aU (77,1;) e L2(o,T; Wo"''^(f})) x L'^ ,T -Wf^ {^)) . We make clear that if 
{g,w,u) is a triple satisfying (|2.ip and (14. ip . then the pair {g,u) is also a weak 
solution according to Definition 12. II 

To obtain a stable numerical method, the mixed finite element formulation of 
(|4.ip is posed with the velocity Vh in a div-conforming space 
and vorticity Wh in a curl-conforming space Wh{n) C Wo"''^(0). There exists 
several such spaces, however here we will use the Nedelec spaces of first order 
and first kind. We choose these spaces for their simplicity and since the most 
natural choice of approximation space for the density is then the space of piecewise 
constants. We will continue to denote this space by Qh{^)- 

This choice of finite element spaces is also very convenient since they can be 
related through the exact de Rham sequence 

„ C rr gi'^d Tir curl Tr div „ _ 

> Sh > Wh > Vh > Qh > 0. 

Thus, we can use spaces orthogonal to the range of the previous operator, i.e., 

Wl:^ {wh e Whicnvlwh = 0}^ nWh, := W G F,,;divi;,, = Oj^nV^, 

to deduce the decompositions 

Wh = DSh + Wl-'^, Vh = curl Wh + , 

together with the discrete Poincare inequalities 

hhWL^^n) < C'l|divVh|li2(o) , 1 1 -^/i II £2(0) < C' l|curlw^|1^2(o) , 

Consequently, as with the previous method, the mixed finite element method also 
admits Hodge decompositions, which in turn implies that a discrete equation for 
the effective viscous flux can be derived. 

We need the following compactness property of the space Vj^'^ . The proof is 
given in [4', Appendix A] . 

Lemma 4.1. Let {t'/i}/i>o be a sequence in V^^'^ such that \\ divt;^||i2(o) < C, 
where the constant C > is independent of h. Then, for any ^ S K^, 

\\vhix) - Vh{x - e)lli=(o) < + ICn^ll divvhllmn), 

where the constant C > is independent of both h and ^. 
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Definition 4.2 (Mixed finite element method). Let and /™ be as given in 
Definition 1^21 Determine functions 



{Q'f:,wl\ul') e Quin) X Wnin) X Vh{n), m=l,...,M, 
such that for all (ph G Qh{^), 
(4.2) 

/ 4[g^]<l>n dx-AtY, I . v)+ + g™« • ^Y) [<l>h]v dS{x) = 0, 



and for all {'qh.Vh) e Wh{^) x Vh{n), 



/ ^ curl w]^Vh + + A) div it™ - _p(£i")] div Vh dx 
Jn 

/ w^rjh — w™ curlry^ dx = 0. 
Jo 

In (|4.2p . (m;i • i^)+ = max{M;i ■ i^, 0} and (m/j • i^)^ = min{M;i • i^, 0} 



(4.3) 



The existence of a solution to the discrete equations (|4.2p - (|4.3p is proved in 
[4]. Moreover, for any ni — 1, . . . , M, 



p lit IIL 

/ P{e^)dx + Y.^^\\ '"'fillw™rl,2(0) 

fe = l fc = l 

+ E / ^"(^'tt)I^'r'F dx + Y^Y.^^ P"{g'^Mrr\nt->y\ dx 
< / P{go) dx + CYMlftWhin). 

•'^ k=l 



where e [p^ \ p^] and e [p^, g'l]. 

For each fixed /i > 0, the numerical solution {(p™, "if™, 1*")}^^^ is extended 
to the whole of (0, T) x by setting 

(4.4) {gh,Wh,Uh){t) ^ {g';;',w',^,u'^), te(i„_i,i„), m = l,...,M. 

In addition, we set gh{0) = ^?°• The main result in [4] is that the sequence 



{gh, Wh, Uh}h>a converges to a weak solution in the sense of Definition 12.11 

Theorem 4.3. Suppose f e i^((0,r) x Q), go e L^{n) if j > 1, and 
Polog^?o G L^{fl) if J — 1. Let {(p^, lO/j, m/j)}^^^ be a sequence of numerical 
solutions constructed according to (14. 4p and Definition \4-S\ Then, passing if nec- 
essary to a subsequence as /i — > 0, Wh w in L2(0,T;Wo™'^'''(f])), 

ghUh gu in the sense of distributions on (0, T) x Vl, and 
gh ^ g a.e. in (0,r) x fl, where the limit triplet {g,w,u) satisfies the mixed form 
(|4.ip . and consequently {g,u) is a weak solution as stated in Definition \2.1\ 

Comments on the proof of Theorem \4.S\ The proof of convergence to a weak 
solution of the continuity equation (|l.ip is similar to the corresponding step in 
the proof of Theorem 13.31 The difference is that the compactness of the velocity 
approximation now requires a different argument. In particular. Lemma l4.ll must 
be employed. The proof of convergence to a weak solution of the velocity equation 
(|1.2p is also similar to the proof of Theorem 13.31 However, a difference is that 
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the weak sequential continuity of the effective viscous flux can now be obtained by 
using test functions in V^°'^(f7) satisfying divvh = Qh- Strong convergence of the 
density is then obtained as in the proof of Theorem 13.31 

5. Extension to the Stokes approximation equations. 

In this final section we present an extension of the previous finite element 
method to the following system 

(5.1) dtg + diY{gu) ^Q, in (0, T) x 1] 

(5.2) -gdtu- ^i/\u~ XDdiYU + Dp{q) ^Q, in(0,r)xO, 

where P = jnj /f2 denotes the average initial density. The equations (|5.ip - (|5.2p 
is known in the literature as the Stokes approximation equations. The system is 
almost identical to the compressible Stokes system (ll.ll) - (|1.2p . the difference being 
the inclusion of the time derivative term in (|5.2p . 

The finite element method is similar to the mixed method of Section |4l and 
as such it is only applicable to the case of Navier-slip boundary conditions (|1.5p . 
Furthermore, for technical reasons, convergence is proved only in the case "f > 
The method is constructed and analyzed in [6]. 



Definition 5.1 (Numerical scheme). Let p° be as given in Definition 13.21 
Determine functions 

{d^.w^. K) G X X Vhm, m = i, . . . , m, 

such that for all (t)h G Qh{^), 
(5.3) 

/ dx-AtJ2 f ■ + 0+« ■ IMr dS{x) = 0, 

and for all {T]h,Vh) £ Wh{^) x Vh{n), 

/ dt [u]l'']vh + n curl w]^Vh + [{fJ- + A) div - p{gT)] div Vh dx = 0, 

(5.4) 



In 

wl^rjh ~ M™ cmlrjh dx = 0, 
for m — 1, . . . , M. 



Existence of a numerical solution and various properties of these solutions hold 
as in the previous section with only minor modifications. We extend {gjj, w^, Uh}h>o 
for fc = 1, . . . , M to functions {{Qh,wii, Uh)}h>o defined on all of (0, T) x fl as in 
(|4.4p . The main result in [6] is that the sequence {{Qh,Wh,Uh)}h>o converges to a 
weak solution of the Stokes approximation equations (|5.ip - (l5.2p . The notion of a 
weak solution is similar to that in Definition 12.11 



Theorem 5.2. Suppose 7 > ^ and go G L'^{Q) . Let {{gh,Wh,Uh)}f^yQ be 
a sequence of numerical solutions constructed according to Definition \5.1\ Then, 
passing if necessary to a subsequence as ft, — > 0, w in L2(0,T;Wo'"''2(r!)), 

Uh u in L^(0, T; W(5''^'^(r2)), ghUh gu in the sense of distributions, and 
gh — > g a.e. in (0,T) x fl, where the limit triplet {g,w,u) is a weak solution to the 
Stokes approximation equations (|5.ip - (|5.2p . 
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Comments to the proof of Theorem \5.2l The proof of convergence is similar 
to the proof of Theorem 13.31 However, the proof of higher integrabihty on the 
density and the proof of weak sequential continuity of the effective viscous flux now 
requires additional arguments in order to handle the time derivative in (|5.4p . In 
particular, the continuity scheme (|5.3p needs to be used to handle the term uvt, 
where Vt G V°'^(r2) satisfies divvt = Qt- For technical reasons, we must then 
require 7 > Moreover, the higher integrabihty estimate for the density now 
gives p{Qh)Qh Gb L^{^,T] L^iyi)). In addition, in this case we need in fact strong 
convergence of the velocity, Uh — > u. This is obtained through an Aubin-Lions 
argument using the spatial compactness on the velocity together with weak control 
of d!l[uh\. Due to space limitations we refer the reader to [6] for details. 

References 

[1] S. Brenner, J. Cui, L-Y. Sung, A nonconforming finite element method for a two-dimensional 
curl-curl and grad-div problem, Numer. Math. 109 4 (2008), 509-533. 

[2] M. Crouzeix, P-A. Raviart, Conforming and nonconforming finite element methods for solving 
the stationary Stokes equations I., RAIRO Anal. Numr. 7 (1973), 33—75. 

[3] E. Feireisl, Dynamics of compressible viscous flow, Oxford Lecture Series in Mathematics and 
its Apphcations, 26 (2004). 

[4] K. H. Karlsen, T. K. Karper, A convergent mixed finite element method for a semi-stationary 
barotropic compressible Stokes system. Preprint. 2008. 

[5] K. H. Karlsen, T. K. Karper, A convergent non-confirming finite element method for a semi- 
stationary barotropic compressible Stokes system. Preprint. 2008. 

[6] K. H. Karlsen, T. K. Karper, A convergent finite element method for the Stokes approximation 
equations. Preprint. 2008. 

[7] P-L. Lions, Mathematical topic in fluid mechanics. Vol 2: Compressible models, Oxford Lec- 
ture Series in Mathematics and its Applications, 10 (1998). 

[8] E. Weinan, Dynamics of vortex liquids in Ginzburg-Landau theories with applications to 
superconductivity, Phys. Rev. B. 50 2 (1994), 1126-1135. 

[9] J. Zhao, D. Hoff, A convergent finite-difference scheme for the Navier-Stokes equations of 
one- dimensional, nonisentropic, compressible flow, SIAM J. Numer. Anal. 31 5 (1994), 1289- 
1311. 

[10] J. Zhao, D. Hoff, Convergence and error bound analysis of a finite- difference scheme for the 
one-dimensional Naviei — Stokes equations, AMS/IP Stud. Adv. Math, 3 (1997). 

Center of Mathematics for Applications (CMA), University of Oslo, P.O. Box 1053, 
Blindern, N-0316 Oslo, Norway 

E-mail address: keiinethk@math.uio.no, t .k.karperOcma.uio .no 



